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We study the statistical properties of time distribution of seimicity in California by means of a new 
method of analysis, the Diffusion Entropy. We find that the distribution of time intervals between 
a large earthquake (the main shock of a given seismic sequence) and the next one does not obey 
Poisson statistics, as assumed by the current models. We prove that this distribution is an inverse 
power law with an exponent fj, = 2.06 ± 0.0L We propose the Long-Range model, reproducing the 
main properties of the diffusion entropy and describing the seismic triggering mechanisms induced 
by large earthquakes. 
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The search for correlation in the space-time distribu- 
tion of earthquakes is a major goal in geophysics. At 
the short-time and the short-space scale the existence 
of correlation is well established. Recent geophysical 
observations indicate that main fracture episodes can 
trigger long-range as well as short-range seismic effects 
[H 13 0- OJ ■ However, a clear evidence in support of these 
geophysical indications has not yet been provided. This 
is probably the reason why one of the models adopted to 
describe the time distribution of earthquakes is still the 
Generalized Poisson (GP) model [a, la> tl ISl IS| • Basically 
the GP model assumes that the earthquakes are grouped 
into temporal clusters of events and these clusters are 
uncorrelated: in fact the clusters are distributed at ran- 
dom in time and therefore the time intervals between 
one cluster and the next one follow a Poisson distribu- 
tion. On the other hand, the intra-cluster earthquakes 
are correlated in time as it is expressed by the Omori's 
law |lOlll]| . an empirical law stating that the main shock, 
i.e. the highest magnitude earthquake of the cluster, oc- 
curring at time to is followed by a swarm of correlated 
earthquakes (after shocks) whose number (or frequency) 
n(t) decays in time as a power law, n(t) oc (t — to)~ p , 
with the ex pon ent p being very close to 1. The Omori's 
law implies |12( that the distribution of the time intervals 
between one earthquake and the next, denoted by r, is a 
power law i/j(t) oc t~ p . This property has been recently 
studied by the authors of Ref. [I^] by means of a unified 
scaling law for V'LjmCt), the probability of having a time 
interval r between two seismic events with a magnitude 



larger than M and occurring within a spatial distance L. 
This has the effect of taking into account also space and 
extending the correlation within a finite time range t* , 
beyond which the authors of Ref. ^1 recover Poisson 
statistics. 

In this letter, we provide evidence of inter- clusters cor- 
relation by studing a catalog of seismic events in Cali- 
fornia with a new technique of analysis called Diffusion 
Entropy (DE) |L3l Il4| . This technique, scarcely sensitive 
to predictable events such as the Omori cascade of after- 
shocks, is instead very sensitive when the deviation from 
Poisson statistics generates Levy diffusion EH • This 
deviation, on the other hand, implies that the geophysical 
process generating clusters has some memory. In Fig. Q] 
we report the sketch of the typical earthquakes frequency 
vs time in the catalog considered. By = U + i — U we 
indicate the time interval between an earthquake and the 
next. Each peak of frequency (cluster) in figure includes 
the time location of a main shock. The time interval be- 
tween one peak and the next is reported in figure and 
is denoted by the symbol r| m ' , where the superscript m 
stands for main shock. We assume that two distinct times 
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rl m l are not correlated, i.e. (r| J rj J ) = ([t^) }5 itj . 

This assumption does not conflict with the departure of 
the distribution of the times from Poisson. On the 
other hand the intercurrence times r^'s are correlated. In 
fact, after a main shock the earthquake frequency is high. 
Consequently, the r's are short and a short value of r is 
followed with a large probability by another short value. 
For the same reason we expect that, far from a main 
shock, and prior to the next one, a long value of r is fol- 
lowed by another long value of r. This implies that the 
correlation function (TjTj) does not decay to zero after 
one step, and that it survives for the whole time interval 
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between two consecutive main shocks. This means that 
the long-time relaxation of the correlation function (t^t.,) 
is determined by the -statistics 0], being faster or 
slower, according to whether i/)(rt m l) is Poisson or not. 
As we shall see, the non-Poissonian condition is straight- 
forwardly assessed by the DE, if ^(r^) has an infinite 
second moment. It is important to stress that the model 
of Fig. 1 will be used to support the results of the paper 
with the study of artificial sequences, but in no way it 
implies an a priori identification of the main shocks for 
the DE method to work. 

The DE method, as almost all the techniques used to 
detect correlation in a time series £(t), is based on the dif- 
fusion process of the auxiliary ir-space through the equa- 
tion: 



dx . , . 



(1) 



In the case under study here, the stochastic variable £(t) 
is constructed by setting £(f) = 1 (or £(t) = M, with M 
being the earthquake magnitude) if an earthquake occurs 
at time t, and setting £(t) = otherwise. In practice £(t) 
is a string of long patches of O's occasionally interrupted 
by l's. With this prescription we build up a diffusion 
process in the x-space |l6j . We construct many distinct 
trajectories, labeled by the integer index n = 1,2, ., ac- 
cording to the prescription 



x n (t) 
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£(t')dt', 
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where the generic walker x n (t) takes a step ahead, of 
either length 1 or M, every time an earthquake occurs. 
Note that At is the resolution time, set in our case to 
be equal to 1 min. All the trajectories x n (t) occupy at 
time t = the position x — and then spread up over 
the x-axis as a result of their partial or total random 
nature. We study the probability distribution p(x,t)dx 
of finding the walker position at time t in the interval 
[x, x + dx]. In the reference frame moving with velocity 
v = W, where W is the average number of earthquakes 
in the time interval At, the diffusion process is expected 
to fulfill the scaling condition: 



p(x,t) 
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where F(y) is a function with a form dictated by the 
statistics of the process, and 6 is the so-called scaling 
parameter. According to the results of Ref.[l3l Il7| the 
evaluation of the scaling parameter S requires the use of 



S(t) = — J_ dx p(x,t)\n[p(x,t)]. In fact, using Eq. 1(3), 
we get after some simple algebra: 



S(t) = A + 5 ln(t). 



(4) 



This means that the entropy of the diffusion process is 
a linear function of ln(i) and a measure of the slope is 
equivalent to the determination of the scaling parameter 
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FIG. 1: A sketch of the typical earthquakes frequency vs 
time. In correspondence to each main shock we observe a fre- 
quency peak determined by the after shock swarm. The peaks 
decay according to the Omori's law. The horizontal arrows 
indicate the time intervals r| m ' between two consecutive main 
shocks. The DE method gives information on the distribution 
of these time intervals. 



r5. It is now important to observe that the DE method 
has the interesting property of detecting the statistics of 
really random events, as recently discovered by the au- 
thors of Ref. 0,0]. We refer the reader to these papers 
for mathematical proofs. For the purposes of this letter, 
the following remarks should suffice. The time intervals 
t's are correlated, as shown with the help of Fig. ^ 
while the DE method rests on the Shannon entropy and 
the Shannon entropy increases only as a consequence of 
the occurrence of really random events. For the time be- 
ing we rule out the possibility that the Shannon entropy 
increase is determined by a deterministic bias \Vt^, and 
so by a non-stationary condition. In the stationary con- 
dition, the only source of entropy increase is given by 
the occurrence of clusters of seismic events because the 
T [m]' s are n0 £ correlated. The deviation from the Pois- 
son statistics are easily detected by the DE method if 
t/>(tM) produces anomalous diffusion. Let us consider 
the non-Poissonian waiting time distribution: 



( T [m]) A 



(5) 



The condition /i > 3 implies a finite second moment, and 
consequently ordinary diffusion with S = 0.5. The con- 
dition 2 < /i < 3, on the contrary, produces an infinite 
second moment, and, conse quen tly, through the gener- 
alized central limit theorem [l4(, the anomalous scaling 



(6) 
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The condition 1 < ji < 2 produces an anomalous scaling 
with 



S = fi-1 



(7) 
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FIG. 2: The Shannon entropy S(t) of the diffusion process 
as a function of time, in a logarithmic time scale. From top 
to bottom, the curve refers to all events (full circles) and 
to events with threshold M — 2,3,4 (open symbols). The 
straight lines have the slope S = 0.94. 
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FIG. 3: The Shannon entropy S(t) of the diffusion process as 
a function of time, in a logarithmic time scale. Open squares 
and full circles are, respectively, the results of the GP and of 
the LR model. The two straight lines, have the slopes 5 — 0.5 
and 5 = 0.94. 



and would imply non-stationarity (as in the presence of a 
deterministic bias). As we shall see, the non-Poissonian 
statistics of the distance between two clusters (i.e. be- 
tween large events) is detected by the DE method, yield- 
ing the anomalous scaling parameter S = 0.94 ±0.01. 

The catalog we have studied covers the period 1976- 
2002 in the region of Southern California spanning 20° 
N -45° N latitude and 100° W 125° W longitude [H . 
This region is crossed by the most seismogenetic part of 
the San Andrea fault, which accommodates by displace- 
ment the primarily strike-slip motion between the North 
America and the Pacific plates, producing velocities up 
to 47mm /yr j^. The total number of recorded earth- 
quakes in the catalog is 383687 and includes the June 
28 1992 Landers earthquakes (M = 7.3), the January 17 
1994 Northrifge earhquake (M = 6.7), and the October 
16 1999 Hector Mine earthquake (M = 7.1). Geophysical 
observations point out that these large earthquakes have 
triggered a widespread increase of seismic events at re- 
mote distances in space and in time 0, Q . The coupling 
of the sources of stress change (i.e. large earthquakes oc- 
currence) and seismicity triggering mechanisms is a pri- 
mary target of geophysical investigations, and, as shown 
below, is revealed by the DE analysis. 

In Fig. |21 we report the results of the DE method. 
The analysis was performed by setting £(i) = 1 when an 
earthquake occurs at time t (independently whether it is 
a main or an after shock), and £(t) = if no earthquake 
happens. In full circles we plot the entropy S(t) as a func- 
tion of time when all the seismic events of the catalog are 
considered (independently of their magnitude M) . After 
a short transient, the function S(t) is characterized by a 
linear dependence on Int. A fit in the linear region gives 
a value of the scaling parameter S — 0.94 ± 0.01 at 95% 
of confidence level. We next consider (open symbols in 
Fig. |2l only the earthquakes with magnitude larger than 
a fixed value M = 2,3, 4. We see that, regardless of the 
value of the threshold M adopted, the function S(t) is 



characterized by the same long-time behavior with the 
same slope. This indicates that we are observing a prop- 
erty of the time location of large earthquakes. This leads 
us to conclude that the time intervals between two large 
events fit the distribution of Eq. JSJ) , with the value of /x 
related to 5 through Eq.®, = 2.06 ± 0.01. In fact 
numerical checks on the time series under study have 
supported the stationary assum ptio n and ruled out the 
alternative condition of Eq.lQ) 22] . Our conclusion is 
also supported by the use of two different walking pre- 
scriptions. The former rests on assuming £(i) equal to 
the magnitude M of the earthquake, at each time when 
an earthquake occurs. The latter sets with equal prob- 
ability either £(t) = 1 or ((t) = —1 when an earth- 
quake occurs . Both methods yield the same exponent 
fx = 2.06 ± 0.01. 

The results of our statistical analysis supports the geo- 
physical arguments that earthquakes of large magnitude 
produce strain diffusion. Unlike co-seismic deformation, 
which is practically instantaneous, the strain diffusion 
ensuing an earthquake of large magnitude produces post- 
seismic stress changes, generating a remarkable increase 
of rate of seismicity at locations hundred of kilometers 
away and over time span up to several years. Consis- 
tently the distribution of the time intervals between two 
large earthquakes is not a Poisson function. Our con- 
clusion is reached under the important assumption that 
the sequence {t- } is not affected by any deterministic 
bias, this being a possible source of ballistic scaling [T^ |. 
If this condition applies clusters occurrence would show 
deterministic trends on the scale of the whole sequence. 
This is an attractive possibility that does not seem to be 
ruled out by the current literature on this subject [23| . 
However, here, we adopt the explanation that the wait- 
ing time distribution is given by Eq. I® with [i > 2, this 
being the unique consequence of the stationary assump- 
tion. 

We now illustrate how the DE method works on two 
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artificial earthquakes time series: the first generated by 
means of the GP model, and the second generated by a 
new model, the Long-Range (LR) model, that we pro- 
pose as a better model to reproduce the properties of the 
catolog considered. In the LR model the earthquakes 
are grouped into temporal clusters, and, as in the GP 
model the number of earthquakes in a cluster follows the 
Pareto law, i.e. a power law distribution with exponent 
equal to 2.5 HHG|. 

The events within the same cluster 
are distributed according to the Omori's law: the inter- 
val t follows a power law with exponent p = 1. How- 
ever, in the LR model the time distance between 
one cluster and the next follows a power law with ex- 
ponent /i = 2.06, rather than a Poisson prescription as 
in the GP. Notice that this value of fi is close to the 
border between stationary and non-stationary condition 
[T^ . The two sequences have the same time length. We 
choose the number of clusters in order to have the same 
total number of earthquakes as in the real data [24| . The 
result of the DE on the artificial sequences is reported in 
Fig. |21 The GP model is characterized by a long-time 
behavior that, as expected, fits very well the prescription 
of ordinary statistical mechanics, with 5 = 0.5. The LR 



model yields the quite different scaling 5 = 0.94. It is 
also clear that the LR model yields a behavior qualita- 
tively similar to that produced by the real data of Fig. 
2 as well as the same scaling parameter 6 = 0.94, while 
the GP fail reproducing both properties. 

In conclusion, this paper is the first application of the 
DE method to study the statistical properties of earth- 
quakes time distribution. We have found that there ex- 
ists a correlation mechanism beyond the Omori's law. 
Both intra-cluster swarms and inter-cluster distances 
obey an inverse power law prescription, the former be- 
ing i/>(t) cx t- 1 and the latter V(r [m] ) oc ( T ^)-^ with 
fi = 2.06 ± 0.01. We have proposed a new model, the 
LR model, better than the GP model in reproducing real 
data. The method proposed is based on the fact that the 
asymptotic properties of diffusion process generated by 
the seismic events are scarsely sensitive to the memory 
stemming from the Omori's law. They are, on the con- 
trary, sensitive to the anomalous statistics generated by 
the non-Poissonian nature of the time distance between 
two consecutive large earthquakes. 

PG acknowledges support from ARO, through Grant 
DAAD19-02-0037. 
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